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ABSTRACT 

Baryon acoustic oscillations (BAO) provide a robust standard ruler with which to measure the 
acceleration of the Universe. The BAO feature has so far been detected in optical galaxy surveys. 
Intensity mapping of neutral hydrogen emission with a ground-based radio telescope provides another 
promising window for measuring BAO at redshifts of order unity for relatively low cost. While the 
cylindrical radio telescope (CRT) proposed for these measurements will have excellent redshift reso- 
lution, it will suffer from poor angular resolution (arcminutes at best). We investigate the effect of 
angular resolution on the standard ruler test with BAO, using the Dark Energy Task Force Figure of 
Merit as a benchmark. Wc then extend the analysis to include variations in the parameters character- 
izing the telescope and the underlying physics. Finally, wc optimize the survey parameters (holding 
total cost fixed) and present an example of a CRT BAO survey that is competitive with Stage III 
dark energy experiments. The tools developed here form the backbone of a publicly available code 
that can be used to obtain estimates of cost and Figure of Merit for any set of survey parameters. 

Subject headings: cosmology — large scale structure of universe — baryon acoustic oscillations - 
standard ruler test — 21cm intensity mapping 



1. INTRODUCTION 

A standard ruler test with Baryon acoustic oscillations 
(BAO) is considered the most robust and systematics- 
free method to probe the dark energy equation of state 
(jAlbrecht et al.l 120061) . The sound waves which prop- 
agated through a mixture of photons and baryons in 
the early Universe left a distinct oscillatory signature 
in the Cosmic Microwave Background (CMB) and the 
large scale structure of matter (or galaxies) with a char- 
acteristic scale set by the sound hori z on at the epoch 
of recombination ([Peebles fc Yul 119701: iHoltzma r] [l98l 
IHu fc Sugivamalll996l; lEisenstein fc Hulll998ft . Measure- 
mcnts of the CMB provide this sound horizon scale, 
and we can then use the BAO in low-redshift cluster- 
ing as a standard ruler to measure the distances to var- 
ious redshifts , and therefore the equation of state of 
dark energy (IHu fe White! 119961 lEisenstein et al.l 11998 



Eisensteinll2003l: iBlake fc Glazebrookll2003H Linderl 
Hu fc Haimanll2003l : iSeo fc Eisensteinll2003| ) 
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While detections of BAO so far have used the opti- 
cal band of the electromagnetic spectrum, using spectro- 
scopic or photometric techniques, there has been grow- 
ing interest in the feasibility of using 21cm emission from 
neutral hydrogen (IWvithe fc Loebll2008l: lLoeb fc Wvithd 



120081 IWvithe fc Loebll2009HVisbal et ai.ll200"9j) to d etect 
BAO at high dMao et all 120081: IWvithe et al.ll2008fl and 
medium (jChang et al.ll2008HAnsari et al.ll2008f ) redshifts. 
Observing BAO from intensity mapping of 21cm emis- 
sion, especially at medium redshift near z ~ 1, has sev- 
eral advantageous features. The BAO is a relatively weak 
feature on large scales, and therefore very large survey 
volumes are needed. A ground-based cylindrical radio 
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telescope (CRT) with a large field of view can easily cover 
most of the sky in a relatively short time. Second, the 
electronics required for frequencies near v ~ 1GHz are 
cheap and easy to build. Digital electronics with precise 
timing offer high precision (better than ppm) frequency 
and, hence, redshift measurements. Third, radio surveys 
rely on different tracers of large scale structure (neutral 
hydrogen) than do optical survey (luminous galaxies). 
Seeing the signal in two different sets of mass tracers 
would be compelling evidence for its robustness. 

A disadvantage of radio waves is their longer wave- 
lengths, making it more difficult to obtain high angular 
resolution. For example, a 100m by 100m telescope will 
have an angular resolution of (A/L) ~ 15 arcminutes at 
z ~ 1. Radio observations of cosmological hydrogen are 
severely contaminated by foreground sources, primarily 
Galactic synchrotron radiation. We will not consider the 
problem of foreground removal here, but will show that 
foregrounds are smooth at the BAO scale and then as- 
sume that they can be eliminated to the accuracy of the 
statistical errors. 

In this paper, we consider a ground-based 21cm in- 
tensity mapping with a CRT targe ting a redshif t rang e 
of 0.2-2, motivated by the study of iChang et al.l (|2008l ). 
In § [3J we derive the signal-to-noise of the power spec- 
trum for a given set of radio telescope configuration pa- 
rameters, discuss galactic shot noise and the effect of 
foregrounds, and explain the assumptions used in our 
method of Fisher matrix projections. We model the ef- 
fects of angular resolution with a window function that 
approximates the effects of a typical CRT telescope re- 
sponse function on the power spectrum. Using both 
Fisher matrix analysis and Monte Carlos, we study the 
issue of angular resolution in detail in §|3l focusing on the 
impact of resolution on the determination of dark energy 
parameters. Then in § HI wc consider a simple compact 
CRT telescope, and investigate the Dark Energy Task 
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Force (DETF) Figure of Merit (FoM) (|Albrecht et al.l 
1200(1 as a function of telescope parameters. The main 
result of the paper is captured in Fig. [6] where the FoM 
is plotted as a function of these parameters. In § [5j we 
optimize the telescope configuration parameters for the 
maximum FoM and show that an inexpensive CRT BAO 
survey can be competitive with Stage III dark energy 
exp eriments (for optimazation st udies in optical surveys, 
see iParkinson et al.l 120071 l2009j) . In § [6] we summarize 
the results. 

2. INGREDIENTS FOR PROJECTIONS IN A 
RADIO TELESCOPE 

In this section, we derive the signal-to-noise ratio in 
the power spectrum for a given configuration of a cylin- 
drical radio telescope. Once the signal-to-noise ratio per 
Fourier mode is determined, it is straightforward to prop- 
agate this to the errors on the acoustic scale and ulti- 
mately on dark energy parameters. 

2.1. Signal to Noise Ratio 

In a radio telescope, the measurement is power (i.e., 
Watts) received by the antenna, and we denote this with 
a lower case p to distinguish it from the power spec- 
trum P(k) in wavenumber space. The average and the 
standard deviation of power received by an antenna per 
bandwidth Af due to the instrumental and sky noise are 



Pn = k B (gT sky + T a )A/, 



-Pn, 



(1) 
(2) 



where k B is Boltzmann's constant, g is the gain, T sky is 
the average sky temperature (e.g., due to foregrounds), 
and T a is the average antenna noise temperature or the 
amplifier noise temperature. After M such measure- 
ments, the uncertainty associated with p^ per pixel for 
a compact array telescope is 



k B (gf sky + f a )Af = k B (gf sky + f a )Af 

VU^AJ 



, (3) 



where ij nt is the integration time per pixel per band- 
width. 

Meanwhile, the BAO signal is the spatial tempera- 
ture variation due to the clustering of neutral hydrogen, 
<5 ps (n, z) = psO~Hi{n, z) where Sni(h, z) is the fractional 
overdensity of neutral hydrogen at angular position h and 
at redshift z. The average power (again in, e.g., Watts) 
due to the 21cm line is 



Ps = k B gT sig Af, 



(4) 



with T s i g the average brightness tempe rature due to the 
21cm line. This has been estimated (|Barkana fc Loebl 
[2007t iPritchard fc Loebll20u8t iChang et al.ll2008[) to be 
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xm(z)fl Hfi h(l + zf 
H(z)/H 



mK, 



(5) 



where xhi is the neutral hydrogen fraction at z, ^ln,o 
is the ratio of the hydrogen mass density to the critical 
density at z = 0, and H(z) is the Hubble parameter at 
z. 



The measured power spectrum of the 21cm intensity 
P(fc,/i) will include both signal and noise: 



P(k, / i)= P 2 s P m (k,ti) + V R at 



(6) 



where Phi(A:,^) is the power (in units of /i~ 3 Mpc 3 ) due 
to the large scale structure of the neutral hydrogen (i.e., 
the signal) and Vr is the volume of a pixel (also with 
units of /i _3 Mpc 3 ). We present a simple derivation of 
this equation in Appendix A. Since we are interested in 
structure on large scales, we assume that Pm(fc,/i) = 
b 2 P(k,p) where P(k,/i) is the underlying matter power 
spect rum with linear theory redshift distortions (jKaiserl 
\x is the cosine of the angle between the wavevector 
k and the line of sight, and b is a constant bias factor. 

Additionally, there will be a galactic shot noise con- 
tribution due to the discreteness of the HI sources with 
an effective number density n. Adding this noise leads 
to our final expression for the signal to noise per Fourier 
mode: 
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V R 



is the sensitivity per pixel. As will 
we approximate P s hot = 1/n = 



where %±%M 
be explained in § 

100/i~ 3 Mpc 3 . Assuming our fiducial CRT configuration 
as explained in § OH P H i(fc, fi = 0) ~ 900h~ 3 Mpc 3 and 

(gT Bk y+T ) y _ 1800/l -3 Mpc 3 at z „ L 

gJsigVtint A/ J 

For a power spectrum averaged over a wavenumber 
range of width dk, the signal to noise is increased by 
the square root of the number of modes, so 
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where V sur is the total volume of the survey. 

We modify this equation further by including a window 
function W that approximates the effect of an instrument 
response function. 
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(9) 



We discuss details of W in S [23 



2.2. Telescope parameters for our fiducial compact 
array CRT 

The input parameters required to achieve a given signal 
to noise (Eq. |9]) depend on the configuration of the tele- 
scope. The important parameters are the length of the 
cylinder L cy \, the width of the cylinder W cy i, the spac- 
ing of receivers dp (or the number of receivers Nf cc d = 
Lcyi/dp), and the total survey time N yelir . We consider 
a rectangular compact array so that iV cy i = L cy i/W cy \. 
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NcyrLc,i/W cyl 

Nfeed^yl/df 



Fig. 1. — General configuration of a compact CRT array. L cy i is 
the length of the cylinder, W cy i is the width of the cylinder. We 
assume the cylinder is square and compact so W cy \ times the num- 
ber of cylinders is set equal to L cy \. Also the number of antenna 
feeds along a single cylinder is equal to the length of the cylinder 
divided by dp, the spacing of the receivers. 
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Fig. 2. — Expected errors on the power spectrum for our fiducial 
survey in Table [T] We assume a power spectrum at z = 1 in real 
space, after properly accounting for the nonlinear degradation on 
BAO, and the instrumental noise at z = 1 while assuming the 
entire survey volume. A wavenumber bin width of 0.01/i Mpc -1 is 
assumed. 



Fig. Q] shows the general configuration for the compact 
array. For this configuration, 
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2-K-— = 2ttA 
dp L cy i 



(10) 
cos(6)d9 
(11) 



Integration time per pixel, t- mt ss N yeeiT Df 



1 A 



2tt W c 



cyl 



(12) 

Above and throughout, we have made the small an- 
gle approximation, assumed critical sampling, and ne- 
glected side lobes. We also assume a CRT telescope at 
the equator as our fiducial configuration for simplicity, 



TABLE 1 
Fiducial CRT configuration. 





Low rcdshift 


H 


gh rcdshift 


Parameters 


0.66 < 2 < 1.24 


1.22 < z < 2.11 


Length of Cylinder, L cy \ (m) 


99.8 




142.8 


Feed spacing, dp (m) 


0.39 




0.558 


Width of Cylinder, W cyi (m) 


14.3 




14.3 


Duty factor, Df 


0.5 




0.5 


Nycar (years) 


1.40 




0.87 


ZHlf2H,0 


0.00037 




0.00037 


bias 


1.0 




1.0 


Sky temperature, T^y (K) 


10 




10 


Antenna temperature, T a (K) 


50 




50 


gain, g 


0.8 




0.8 


^ahot 


100.0 




100.0 


Note. — For technical reasons, wc assumed 


that 


the redshift 



range is covered by two distinct configurations of the telescope. 

which, according to equation (fTT|) . affords all steradian 
sky coverage in the limit of half-wavelength feed spac- 
ing resulting in coverage to the horizon. We consider 
the perhaps more likely case of a telescope at a lati- 
tude of ~ 30 degrees in section § The angular res- 
olution and the frequency resolution Af determine the 
pixel size, and therefore Vr, and Asurvey determines 
Vsurvey. Note that the effect of A/, i.e., the resolu- 
tion along the line of sight cancels out, as it appears 
both in the sensitivity per pixel and Vr. We neverthe- 
less assume that the frequency resolution is good enough 
(i.e., better than 1MHz, which can be easily achievable 
in this survey) that we can ignore the instrument re- 
sponse function in this direction (sec § 12. 5j) . Wc fix the 
remaining parameters such as the antenna temperature, 
sky temperature, gain, and the observing duty factor to 
be T a = 50K, T sky = 10K, g = 0.8, and D f = 0.5. We 
make a conservative assumption for the neutral hydro- 
gen fraction: xm(z)n„^ = 0.00037 (IZwaan et alj 120051: 
iProchaska et alJl2005t iRao et al.ll2006t ). 



Table [T] lists the telescope parameters for a fiducial 
CRT survey. For technical reasons 0, the rcdshift range 
will likely be covered by two distinct configurations of the 
telescope: in our fiducial survey, the first configuration 
covers 0.60 < z < 1.24 and the second configuration cov- 
ers 1.22 < z < 2.11. These will be combined to obtain 
estimates of the power spectrum in a series of redshift 
bins of size Az = 0.1. The signal to noise on the power 
spectrum will depend on the redshift bin under consider- 
ation. Figure [2] shows the power spectrum at z = 1 with 
associated errors derived from equation (|9]) for the fidu- 
cial survey. To give a sense of the constraining power of 
21 cm surveys, the error bars in this figure assume that 
the total volume of the survey is concentrated at z = 1. 

2.3. Galactic shot noise estimation 

In 21cm intensity mapping we are only able to probe 
the mean 21cm emission in fairly large volumes of space. 
While intensity mapping has the advantage that it in- 
cludes the emission from all the galaxies, no matter 

3 It is difficult to design an array with large fractional band- 
widths (> 50%). By dividing the design into two configurations, 
the fractional bandwidth of our fiducial CRT falls into a reasonable 
range of about 30%. 



4 



how dim, it has the disadvantage that, due to the lu- 
minosity weighting, the random sampling of the density 
field is dominated by a relatively small number of the 
bright galaxies. This sampling noise is known as galac- 
tic shot noise, and is often calculated under assumption 
that the galaxy position and luminosity is randomly se- 
lected with probability proportional to 1 + bS m where 
b is the bias and S m is the mass overdensity. Under 
these assumptions the shot noise adds a scale indepen- 
dent \/n(z) term to the observed overdensity power. Us- 
ing the Schechter function fit to the HIMF (e utral hy- 
drogen mass function) by IZwaan etaLl (I2005T) we find 

™ = Eio 2$a = O.Ol/i 3 Mpc -3 which we use in our 
calculations without evolution. The evolution of n(z) is 
uncertain but shot noise is not the dominant source of 



2.4. Foregrounds 

At the frequencies relevant for this study 21cm emis- 
sion is far from the dominant source of emission. Syn- 
chrotron (plus a smaller contribution of free-free) emis- 
sion will dominate the 21 cm signal by a factor ~ 10 4 . 
This includes emission from our Galaxy and from extra- 
Galactic sources. However as noted by a number of 
authors dDi Matteo et ail 12002b , [Zaldarriaga et all 



iWang et alj2005HMorales et alll 2006; C hang et al 
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synchrotron (and free-free) emission will have a very 
smooth spectrum unlike the BAO wiggles and one can 
therefore distinguish the two. In fact from the inten- 
sity maps one can extract "modes" which will be largely 
uncontaminated by foregrounds. Furthermore these un- 
contaminatcd modes are the ones which are relevant 
for BAO analysis, and this is why one can calculate 
the Fisher information without taking account of the 
foregrounds. This result is suggested in Figure 3 of 
iChang et al.1 (|2008f ) which shows a significant residual 
foreground contamination for k < 0.04ft, Mpc -1 , i.e. on 
scales larger than the BAO scale. Even so their model for 
foreground residuals is extremely pessimistic and one can 
reduce these significantly more than they have assumed. 
To further justify this statement we give an outline of 
how we have set model-independent lower limits for just 
how smooth the synchrotron spectrum must be and also 
present numeric al results. More specifics can be found in 
iStebbinsI (|2010l ). We decompose the observed intensity 
pattern into modes: 



(13) 



Si(k) = / d h / dv m^{n,v) I v {n) 



For the BAO analysis one will want to use modes which 
are localized in the spatial wavenumber, k, describing 
the spatial pattern of 21cm emission. One will further- 
more want these modes localized in angular wavenumber, 
/, allowing one to distinguish modes whose variation is 
primarily radial or angular. 
The synchrotron intensity is given by 







d -^w { v - 

de e ' 



dc 



(14) 



where e = E% \h x B\, E e is the electron energy, 
B the magnetic field, and ^jf- gives the distribu- 
tion of e along line-of-sight n. W is a function 




Fig. 3. — Plotted is the suppression factor, S(K), which is the 
ratio of the maximum amplitude in spectral wiggles in synchrotron 
intensity in units of the mean intensity. It is a function of the log- 
arithmic wavenumber K. In the inset is plotted the logarithmic 
wavenumber corresponding to the 1st BAO peak as a function of 
redshift. Extrapolati ng to get SUCba q) < lCT 100 . These func- 
tions are computed in[Stcbbins ( 20101 V 



that describes synchrotron emission for given electron 
(jRvbicki fc Lightmanl [1985) . This convolution means 
that small-scale wiggles in the distribution function 
show up as wiggles in the spectrum of If, ync . Since W > 
and is smooth one finds that, relative to the mean, wig- 
gles in I* ync are suppressed relative to wiggles in -rf, 
and the suppression is greater as one looks at higher fre- 
quency wiggles. Furthermore since > one finds 

that the amplitude of wiggles in ^jf can be no bigger 
than the mean. Given these constraints it is impossible 
for the amplitude of wiggles in If, ync to be larger than 
a suppression factor times the mean 7^ ync . At the fre- 
quencies of interest the mean J^ ync corresponding to a 
brightness temperature of a few Kelvin so the physical 
bound on wiggles is the suppression factor times a few 
kelvin. The suppression factor will depend on the physi- 
cal scale of the wiggles and is extremely small for wiggles 
which could masquerade as BAO. 

More specifically if one Fourier transforms the inten- 
sity in In j/: I(K) = f™ x I^ nc e~ iKlnu dp then \I(K)\ < 

S(K)I(0) where S(K) > is the suppression factor 
which is plotted in Figurc|3]for small K. S(K) is given by 
the logarithmic Fourier transform of W. At the large K 
relevant to BAO the S(K) is so small that it is difficult to 
determine numerically, but since the Fourier transform 
of C°° functions, such as W, should fall off exponen- 
tially one can extrapolate the S^A'bao) < 10 -100 where 
^bao is the logarithmic wavenumber corresponding to 
the BAO scale. If one chooses modes which limit the 
amount of small- K contamination one can greatly sup- 
press the amount of possible synchrotron contamination. 
One can apply exactly the same argument to free-free 
emission, which depends on the distribution of gas tem- 
peratures in position or e, and obtain similar levels of 
suppression. In practice suppression factors of a 10 -100 
are not achievable due to the contamination from small 
K, and this level is not really needed since this is much 
larger than many other sources of noise such as the pho- 
ton shot noise. However if one understands the beam 
pattern and frequency response of the telescope, one can 
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expect to achieve suppression factors of < 10 for most 
modes on the BAO scale, which is more than enough 
to make foreground contamination a negligible source of 
noise. 

The CRT will be sensitive not only to total intensity 
but also to linear polarization. While at emission the syn- 
chrotron light is linearly polarized with a smooth spec- 
trum, Faraday rotation can cause small scale wiggles in 
the spectrum of each component of linear polarization 
by the time it reaches the CRT. This does not effect the 
total intensity. To avoid possible contamination by Fara- 
day rotation one should combine the radio signals so as 
to make an intensity map with as small a contamination 
linear polarization as possible. 

In this paper, we assume that the foreground is prop- 
erly subtracted, and consider only the noise (i.e, fluctu- 
ations) associated with it. 

2.5. Fisher matrix calculations and instrument 
response 

Once we have estimates of the noise relative to the 
signal for a given survey and instrumental parameters, 
as in equation (0), we can derive error estimates on the 
angular diameter distance Da{z) and Hubble parameter 
H(z) at various rcdshift bins (of size A z = 0.1) using the 
Fisher matrix formalism presented in iSeo fc Eisensteinl 
(|2007| ) : we consider nonlinear degradation on BAO both 
due to nonlinear structure growth and redshift distor- 
tions at each redshift. We combine this with the Planck 
Fisher matrix from the DETF and derive errors on 
dark energy paramete rs and the Figure of Merit (FoM) 
(jAlbrecht et al.l 120061 ). The FoM is defined as the in- 
verse volume of the 95% confidence ellipse in the space 
of wo and w a . Following the convention of the DETF, 
our parameters include the matter density Q m h 2 , baryon 
density Qbh 2 , dark energy fraction fide, curvature frac- 
tion the spectral tilt n s , and the amplitude A, and 
two dark energy equation-of-statc parameters: wq and 
w a . We choos e our fiducial cosmo logy to be consistent 
with WMAP1 dSpergel et al.1120031). 



Recent results (jEisenstein et al 



20071 : ISeo et al1l200ll 



have shown that when the signal to noise ratio per 
Fourier mode is of order two, we can partially undo 
the nonlinear degradation of the BAO signal. For our 
Fisher projection, we assume that we can conduct 
reconstruction and halve the nonlinear damping scale 
of the BAO only when P H i(fc, (J, = 0)/(P N + P sho t) at 
fc = 0.2h Mpc -1 is larger than two. We hereafter denote 
P ffl (fc, fi = 0)/(Pjv + Pshot) at k = 0.2h Mpc^ 1 as nP . 2 . 



While we do not attempt to include an exact instru- 
ment response function, we adopt a reasonable approx- 
imation that describes the general response function of 
CRT in our Fisher matrix analysis. The CRT cylinders 
are oriented in a North-South line so that the field of 
view has a narrow width in right ascension (because of 
the focussing of the cylinders) but a very broad coverage 
in declination. As a consequence the beams formed scan 
the sky at fixed declination and the largest beam separa- 
tion defines a maximum spatial Fourier component (i.e., 
the Nyquist frequency fcNyq = ^/Resolution) that can be 
measured in declination. Higher modes are aliased onto 
lower modes but the alias effect is greatly suppressed be- 
cause of the non-zero beam width. The rotation of the 



earth results in nearly continuous sampling in right as- 
cension. The CRT therefore has some sensitivity to all 
modes beyond fcNyq but the sensitivity decreases around 
&Nyq according to a window function which arises from 
the non-zero beam size. We model the beam shape as- 
suming a uniform illumination of the cylinder aperture 
and calculate the beams that can be synthesized from the 
cross-correlations between cylinders. In both directions 
angular resolution is therefore given by the overall dimen- 
sion of the array. The resulting expression is cumbersome 
but can be adequately represented for our purposes by 
an exponential damping of a signal with a characteris- 
tic damping scale set by the Nyquist limit fcNyq for a 
given angular resolution. As the damping effectively sup- 
presses signals from the modes beyond the Nyquist limit 
along the right ascension, we refer to the effects in both 
directions as Nyquist limit or Nyquist cutoff from now 
on. Such limitation does not apply to wavemodes along 
the frequency direction (i.e., the line of sight) where the 
21cm survey can achieve an excellent redshift resolution. 
For example, a frequency resolution better than 1MHz 
will give a resolution along the line of sight direction 
better than 8/i _1 Mpc even at z = 2. We therefore apply 
no Nyquist limit along this direction. We therefore in- 
troduce the following anisotropic window function W(k) 
and include in the signal to noise calculation in equation 
©: 



W(k) = exv [-1.5(fc x /fc max ) 2 ] 0(fc n 



k y ) (15) 



where O is a Heaviside step function, fc x is a wavevec- 
tor component in the direction of right ascension, fc y is 
a component in the direction of declination, and fc max is 
set to be the Nyquist frequency fcNyq given the angular 
resolution. In a real survey, an instrument response will 
generate a window function with more complex features 
than described above. However, note that the effect of 
the instrument response on the shape of the power spec- 
trum can be precisely predicted and therefore corrected 
for in a real survey. We use equation [9] and [15] as our 
default for § [4] We explain the effect of including this 
high fc cutoff in more detail in § [3j 

An interferometer that forms beams using only cross- 
correlations also has a low fc cut-off that depends on the 
instantaneous field of view along the drift scan direction 
(i.e., in right ascension) but only on the plane fcD CC (= 
k v ) = 0. We have not explicitly corrected for the reduced 
sensitivity at low fc due to the limited field of view, but 
the resulting effect on the total number of modes per fc 
is minimal for the array geometries we consider, so we 
can safely ignore this effect. We note, in addition, that 
the beams that are formed can generally contain signals 
from more than one direction. In our analysis, note that 
we assume that the antenna pattern cuts off sufficiently 
quickly outside the region where the beams are formed 
so that the signal from side lobes are attenuated to a 
negligible level. 

3. EFFECT OF ANGULAR RESOLUTION 

Unlike optical galaxy surveys where the angular resolu- 
tion is of order arcseconds, the expected angular resolu- 
tion of a CRT is a few arcminutes at best. Such a limited 
angular resolution might potentially impair the acoustic 
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Fig. 4. — The effect of introducing the Nyquist frequency limit fc max (i.e., W in eq. [9] and 1151 ) into the Fisher matrix calculation to 
account for a limited angular resolution. Left: Fisher matrix estimates of errors on the distance scale as a function of fc max for different 
cases of nPo.2 assuming a BAO feature at z = 1 in rcdshift space. The arrows point at fcmax = &Nyq for the resolution of 16/i — 1 Mpc 
and 8/i -1 Mpc. Middle: the difference in FoM between including (solid points) and ignoring the Nyquist frequency limit (open points) as 
a function of L cy \, i.e., as a function of an angular resolution for our fiducial CRT configuration. The squares are with the Planck prior 
and triangles are with the Planck + Stage II priors. Right: the corresponding errors on the distance scale Dy at each redshift bin of the 
fiducial CRT. The solid points are with the Nyquist frequency cutoff and the open points are without the cutoff. The same colored points 
in the middle and the right panels correspond to the same value of L cy \. 



peak measurement and therefore distance measurements 
in various ways. The angular resolution of a CRT affects 
the pixel volume Vr in equation (j9)), but also determines 
the range of the wavenumbers available due to the cutoff 
near the Nyquist frequency as modeled with the window 
function W . While the former effect is obvious, the latter 
has not been investigated previously. Before we move to 
the DETF FoM for a CRT, we investigate the effects of 
this window function on the acoustic peak measurement 
and test using Fisher matrix analysis and Monte-Carlo 
simulations. 

We focus on the following effects of angular resolution 
(i.e., the window function). 

• The increase in errors due to the wavenumber cut- 
off near the Nyquist frequency ^Nyq- We can test this 
analytically within the framework of Fisher matrix. 

• Instability of the \ 2 fitting due to the short 
wavenumber range available below k-^ yq . We test this 
using Monte-Carlo simulations. 

First, limited angular resolution will decrease S/N not 
only by increasing Vr but also by limiting the wavenum- 
bers available. Note tha t the Fisher matrix method from 
iSeo fc Eisensteinl (|2007l ) assumes that the power spec- 
trum up to a large enough wavenumber (i.e., fc max = 
0.5ft Mpc -1 ) is available from observation, implicitly re- 
lying on an optical/IR survey. However, the resolution 
of radio surveys will likely limit fcN yq , therefore fc max , 
to be much below 0.5ft Mp c -1 . Due t o the degradation 
of B AO by Silk damping (ISilkl I1968D and nonlinearity 
e.g. iMeiksin. White, fe Peacockfl999HSeo fe Eisensteinl 



200ri Ueong fe Komatsul 120061 lEisenstein et all 120071; 



Crocce fe Scoccimarrol 2008; Matsi ibaral 1200 



0.3 — 0.5ft Mpc -1 contains very little relevant informa- 
tion, but l osing scales k < .3h Mpc -1 may be more 
damaging. iChang et al.l (|2008l ) mention that a Nyquist 
sampled map with a pixel size of ~ 18ft -1 Mpc, and 
therefore fc max = fcN yq = 0.17ft Mpc -1 , is required for 
a BAO survey at z ~ 1, as the information on smaller 



scales is damped away due to nonlinearity. We want to 
probe this quantitatively: how much residual BAO infor- 
mation lies beyond k ~ 0.17ft, Mpc -1 ? This question be- 
comes especially pressing if reconstruction is anticipated, 
in which case we do not want the scales of interest for 
reconstruction to lie beyond the Nyquist frequency. 

We include the effect of the Nyquist frequency limit by 
introducing W as shown in equation © and (TT5|) . which 
effectively limits t he integration range of t he Fisher ma- 
trix calculation in ISeo fc Eisenstein (2007J). We test the 
effect of including this limit as a function of fc max . The 
left panel of Figure 0] shows the Fisher matrix estimates 
of the errors on the the distance measurement Dy as a 
function of fc max for different cases of nPo.2 assuming a 
BAO feature at z = 1 in redshift space; we assume a 
volume of 1ft" 3 Gpc 3 . With fc max ~ 0.2ft, Mpc -1 corre- 
sponding to a resolution of 16ft -1 Mpc (i.e., 23 arcmin at 
z ~ 1), compared to k max > 0.3ft Mpc -1 , the effect of in- 
troducing W is small when the instrumental noise is neg- 
ligible (i.e., nfo.2 3> 10). However, for a realistic survey 
with nPo.2 < 1-5 the error increases by ~ 16% relative to 
an estimate which includes small scales. This degrada- 
tion in errors will be even more severe at high redshift, 
where the linear regime extends to smaller scales leaving 
even more information unobtainable at fixed telescope 
size. 

The survey we consider is composed of a broad range of 
rcdshifts and therefore we observe an integrated effect of 
the Nyquist frequency limit over a range of redshift. The 
middle panel of Figure 0] compares the FoM for our fidu- 
cial CRT configuration as a function of L cy \, i.e., angu- 
lar resolution, with and without considering the Nyquist 
frequency limit, while the rest of telescope parameters 
are held fixed. The right panel shows the effect of the 
Nyquist frequency limit on the errors of Dy at the in- 
dividual redshift bins that comprise the fiducial survey. 
The effect is indeed larger at high rcdshift. Note that 
such an effect would have been much more severe than 
shown here at high redshift if not for the large survey 
area/ volume available at high redshift; if we assume at 
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Fig. 5. — The stability of the \ 2 fitting (i.e., the differences be- 
tween the best fits from different parametrizations) as a function 
of fcmax- The acoustic scale ol is measured from the \^ fitting to 
one Jackknife subsample (i.e., 480 h~ 3 Gpc 3 ) and the differences 
of a in the unit of the true acoustic scale arc shown as a func- 
tion of the parametrization for the broadband shape: for example, 
B2AA means that we used a second order polynomial in k for a 
scale-dependent bias B(k) and a fourth order po lynomial in k for 
a smooth additive term A(k) in Itjeo et al.l H2008T ). Note that the 
fits using fcmax 0.3/i Mpc ^ display robust behaviour regardless 
of the parametrization, but when fc max ~ 0.2/i Mpc- 1 the best fit 
varies with parametrization. 

z = 2 the same volume as at z = 0.7, the errors at z = 2 
will be ~ 1.7 times larger than shown here. Overall the 
effect of the Nyquist limit is small for a CRT array with 
a length larger than 100 ~ 150 meters (i.e., angular res- 
olution better than 10 — 15 arcmin at z ~ 1): <ll%in 
FoM with a Planck prior and even less with an addition 
of the DETF Stage II as priors. A smaller telescope will 
suffer more from the Nyquist frequency limit (magenta 
points). 

The second potential problem with limited resolution is 
related to the way the BAO feature is typically extracted 
from the power spectrum. To make the extraction as 
robust (systematics-free) as possible, the BAO feature 
is extracted by marginalizing over the broadband shape 
of the power spectrum. Poor angular resolution, which 
leads to a limited range of wavenumbcrs available, makes 
the correct estimate of the broad-band shape difficult, 
and this can bias the acoustic scale measurement, and 
therefore distance measurements, as well as increase the 
errors associated with the measurement. 

To test the second problem, we generate 61 random 
Gaussian density fields of 8ft -3 Gpc 3 (i.e., a to tal of 
488ft ~ 3 Gpc 3 ) with 512 3 density grids using a code (jSirkol 
l2005h which produces initial conditions using second or- 
der Lagrangian perturbation theory. To account for non- 
linear degradation of the BAO signal without evolving 
the density field, we smooth the input power spectrum 
based on the expected degradation at z = 1 in real space. 
We Fourier-transform the ensuing distribution, compute 
the power spectrum, and conduct a x 2 analysis to mea- 
sure the acoustic scale using Jackknife resam pling. We 
follow the details of the \ 2 analysis presented in lSeo et al.l 
(|200l . 

In Figure [5l we use one Jackknife subsample (i.e., 480 
ft -3 Gpc 3 ) and show the differences between the mea- 
sured acoustic scales a as a function of the number of fit- 
ting parameters for the broadband shape. For example, 



B2AA means that we used a second order polynomial in k 
for a scale-dependent bias B(k) and a fourth order poly- 
nomial in k for a smooth additive term A(k). A robust 
extraction of a BAO signal is characterized by a best fit 
a that is independent of the polynomial order for a range 
of reasonable choices of polynomials where "reasonable 
choices" means enough flexibility to account for the un- 
known broadband shape while not mimicking the BAO. 
From the figure, while k max > 0.3ft. Mpc -1 shows robust 
behaviour, when /e max ~ 0.2ft, Mpc -1 the best fit varies 
with parametrization: the standard deviation of the best 
fit values among these four parametrizations for the given 
subsample is 7 times larger for /c max = 0.2ft, Mpc -1 than 
for fc max = 0.3ft Mpc -1 . 

We find that the dispersion among many Jackknife sub- 
samples (i.e, the precision associated with the measure- 
ment of the acoustic scale) is also sensitive to a choice 
of parametrization for k mllx < 0.2ft, Mpc -1 . That is, for 
large values of fc max , i.e., the higher resolution data, the 
choice of fitting formula is simple: all reasonable choices 
give the same result. However, in the case of poor reso- 
lution, the marginalizing process over the smooth com- 
ponent of the power spectrum becomes more challenging 
and the fit becomes highly dependent on the fitting for- 
mula. This sensitivity itself speaks to the peril of a low 
resolution survey. 

In summary, angular resolution at the level of 20 ar- 
cmin (i.e., ~ 16ft Mpc) at z = 1 not only increases the 
noise of the survey but also limits the small wavenum- 
ber region where a residual BAO exists. The effect is 
more damaging at high rcdshift where the linear scales 
of interest correspond to smaller angular scales; still the 
impact on the DETF FoM is only 11% for our fiducial 
survey. Second, with a limited wavenumber region as 
a result of the poor angular resolution, it is difficult to 
construct a robust standard ruler test, and the extracted 
acoustic scale will be sensitive to the choice of fitting 
formula used to marginalize over the broadband shape. 
This means that we will need to know the broadband 
shape better, which diminishes the advantages of using 
the BAO as a standard ruler. 

4. FOM AS A FUNCTION OF TELESCOPE 
CONFIGURATIONS 

In this section, we investigate the dependence of the 
final FoM values on various telescope configuration pa- 
rameters. From equation ([9]), it is evident that, at a 
given T s i g , the noise depends on the angular resolution, 
the integration time per pixel, the volume of the survey, 
and therefore on L cy \. N yoar /W cy i, and the rcdshift range 
and dp, respectively. The signal, which scales as T s - lg , will 
vary depending on our assumptions about Qm, and clus- 
tering bias in Phi • In Figure [6j we show the dependence 
of the FoM from a CRT on each of these parameters 
while the rest of the parameters are fixed. We plot two 
values of FoM: first, when we combine the CRT with the 
Planck priors (square points), and second, when we add 
the DETF Stage II survey (cross points). The a;- axis is 
scaled such that the fiducial model parameter presented 
in Table [T] corresponds to unity. The lower panels show 
nP .2 = Phi/(Pn + Phot) at k = 0.2ft Mpc -1 , which in- 
dicates the amount of signal relative to the noise at our 
scale of interest, and <jd v , the resulting error estimate on 
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Fig. 6. — FoM values as a function of various parameters (top panels). The square points are for FoM from CRT+Planck, and the cross 
points are for CRT+Planck+Stage II. The x-axis is scaled such that the fiducial model parameter presented in Table [T] corresponds to 
1 in x-axis. The lower panels show nPo.2 = Phi/{Pn + ^shot) at fc = 0.2h Mpc^ 1 , which indicates the amount of signal relative to the 
noise at our scale of interest, and <Trj v , the resulting error estimate on the distance scale Dy at each redshift bin. The corresponding 
configuration between the top and the bottom panels is denoted with the same colored point. For example, black points correspond to our 
fiducial configuration. 



the isotropic distance scale Dy (Eiscnstcin c t alj|2005l ). 
at each redshift bin for the corresponding configuration 
denoted with the same colored point in the top panels. 
The error on v reflects both V^ ur and nPo.2 as well as 
redshift distortions. 

With our fiducial survey (see Tabled]), we find FoM = 
21 without Stage II, and FoM = 141 including Stage 
II, which is ~ 3 times better than Stage II alone, and 
somewhat worse than Stage II + Stage III. From the 
bottom panels of the figure, one sees that our fiducial 
survey is noise-dominated at high redshift bins (nPo.2 is 
less than one), mainly due to poor angular resolution. 
This results in a larger <jd v in the high redshift bins, an 
effect somewhat offset by a larger volume per redshift bin 
at higher redshift. 

As expected, the FoM increases with increasing L cy \ 
due to increasing (i.e., better) angular resolution. In- 
creasing L cy \ increases nPo.2 and decreases od v over all 
redshift bins. While the improvement on <jd v is in gen- 
eral larger at high redshift, at the lowest redshift bins, 



nPo.2 is now big enough to pass the reconstruction re- 
quirement of rip). 2 > 2. The contribution from these 
lowest redshift bins is very important for improving the 
FoM as L cy \ increases. 

Increasing W cy \ decreases the amount of time that a 
celestial object will sit in the cylinder beam as it drifts 
across the cylinder width, i.e., the integration time, and 
therefore decreases the FoM. Note that decreasing W cy \ 
means increasing the number of cylinders and therefore 
increasing the number of feeds in the case of a compact 
array. Increasing and decreasing 7V yoar will affect the in- 
tegration time in an opposite manner to W cy i, and there- 
fore the two can be adjusted to balance the overall cost 
of the survey and the time required to complete it. 

For the survey we consider here, the receiver spacing 
<If determines the main lobe beam angle, therefore, the 
area of the sky covered. As dp reaches 0.17m, i.e., half 
of our fiducial spacing, the survey area increases to the 
whole sky, and no further improvement in FoM is ob- 
tained, as shown in the right panels of Figure El 
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Fig. 7.— Optimization of the CRT: the maximum DETF FoM 
at given cost. We assume a CRT at a latitude of 35 degrees. 

The figure also shows that the neutral hydrogen frac- 
tion and the clustering bias of the neutral hydrogen are 
important uncertainties in estimating the performance of 
the CRT. Increasing the two parameters has two effects: 
first, the signal increases, therefore increasing the signal 
to noise. Second, some of the low redshift bins now can 
satisfy the nPn.2 > 2 threshold for reconstruction, and 
therefore we can assume a smaller nonlinear degradation 
on the BAO feature for those redshifts. The clustering 
bias of the neu tral hydrogen has been investigated i n 
several paper s dMever et al.l l2007t iWvithe et at! 120091 ) . 
According to iMarin et a l. ( 2009) , large-scale bias of the 
neutral hydrogen is near unity at z = 1, but increases to 
~ 1.3 at z ~ 2. If we include this bias evolution in our 
fiducial model, we find a slightly larger FoM: 24 without 
Stage II and 152 with Stage II for our fiducial survey 
parameters. 

5. MAXIMIZATION OF FOM 

In designing a future CRT array, we want to optimize 
the telescope configuration for the best performance (i.e., 
a largest FoM) while holding the cost fixed. In this sec- 
tion we describe an online calculator which performs this 
optimization. 

For a given survey area, the figure of merit increases 
as the pixel noise and pixel size is reduced. However, 
reducing the pixel noise and pixel size increases the com- 
plexity of the CRT array. For example, while we have 
considered a square compact CRT array in the previ- 
ous sections, we may want to sparsely locate cylinders to 
be more cost-effective, as the cost depends on the num- 
ber of cylinders. The sensitivity per pixel, and therefore 
the noise per pixel, then will also depend on the pack- 
ing factor, i.e., the ratio of the number of cylinders to 
the number of possible locations of the cylinders. We 
can describe a cost function that parametrizes the com- 
plexity of the CRT array. It is not intended that these 
costs include everything that would arise in designing 
and building a large radio telescope, such as site prepa- 
ration, non-recoverable engineering costs, overhead, con- 
tingency etc. These costs should only be used in trying 
to compare sets of design parameters. The cost function 
can be broken into three components; the cost of the dig- 
ital electronics, the cost of the cylinder feed line, and the 
cost of the reflector surface. 

C T = N f N c R e + N f N c d F Rf + N f N c d f W? yl R r (16) 



where Nf is the number of feeds per cylinder, N c is the 
number of cylinders, dp is the feed spacing, and W cy i is 
the width of the cylinders. In addition, R e is the cost per 
digital channel, Rf is the feed line cost per unit length, 
and R r is the reflector cost per unit volume. These three 
cost rates can be scaled from the cost of existing proto- 
types. 

An optimization tool for the maximum FoM per given 
cost that takes into account a large number of CRT engi- 
neering p arameters was developed and is available pub- 
licly at |http : //astro . f nal . gov/21cm| Figure [7] shows 
the output of a design optimization tor various CRT engi- 
neering parameters. Because fractional changes in array 
sensitivity and resolution drop off as the array cost in- 
creases, the corresponding figure of merit levels off as 
well. However, it should be noted that a FoM of 300 can 
be achieved with an array costing $15 million according 
to our cost model. 



6. CONCLUSION 

In this paper we have considered a BAO survey with a 
ground-based 21cm intensity mapping with a CRT tar- 
geting a redshift range of 0.2-2.0. We have tested the 
feasibility of the survey while paying more careful at- 
tention to the effect of limited angular resolution. We 
summarize the results in this paper. 

First, we have tested the effect of angular resolution 
using the Fisher matrix calculations and Monte Carlos. 
The angular resolution at the level of 20 arcmin at z ~ 1 
not only increases the instrumental noise beyond the sig- 
nal but also notably limits the wavenumber range avail- 
able for the standard ruler test. The latter effect further 
increases errors on the acoustic scale, while such effects 
on FoM are small for a CRT array with a size of 100-150 
meters (i.e., an angular resolution of 10 — 14 arcmin at 
z ~ 1). The limited wavenumber range also makes the 
marginalizing process over the smooth component of the 
power spectrum to isolate the BAO more challenging: 
the result highly depends on the fitting formula. There- 
fore, in the case of poor resolution, the data are hard to 
fit reliably. 

Second, we have investigated the dependence of the 
DETF FoM on various telescope parameters while as- 
suming a simple compact CRT array. As expected, the 
FoM strongly depends on the telescope parameters as 
well as the neutral hydrogen fraction and the clustering 
bias. 

Third, we have relaxed the assumption of the compact 
CRT array and have optimized the telescope configura- 
tion for the maximum FoM while holding the cost fixed. 
We find that we can achieve FoM of 300 (with Planck and 
Stage II priors) for a reasonable cost. The FoM/cost cal- 
culator and the optimization tool for a CR T is publicly 
available in |http : / / astro"! f nal . gov/2 £cm[ 
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APPENDIX 

21CM INTENSITY POWER SPECTRUM 

We derive the measured power spectrum in Fourier space, using the following definitions for discrete Fourier trans- 
form: 

% = ^E e *%( f ) ( A1 ) 

X 

5 P {x) = Y,k^ S (A2) 
fc 

^- e «f.(£-/?) =N6jf ff (A3) 

X 

<S^,>V,=P(k)Sf g , (A4) 

where N is the total number of pixels in the map, is the total volume of the map under FFT, and 6 K is the Kronccker 
delta. Here 8 p is the fluctuations in power (i.e., Watts) measured for the pixel due to the large scale structure of HI 
(psShi) and instrument noise such that 

< S p (Xp)5p(Xq) > =Ps(tHl)pq + <4n S W ( A5 ) 

where £hi is the two point correlation function of the neutral hydrogen density field. Parameters ps and ap N are the 
average power due to the 21cm line and the uncertainly associated with the average power p^ due to the instrumental 
and sky noise, respectively, as defined in equation [T] - [4] 
Then, 

< >=< ^E^'^^E 6- ^^) > 



^ E E + «] e**e-*>* 

X q 

A E E E ^PHii^e-^'^-^e^e-^'-^ 



■E p q k 1 ' 



]\[2 E °JW 



2 gi(fc— k'j-x'p 



-a 2 



where Pni(k) is the power spectrum in /i~ 3 Mpc 3 due to the large scale structure of the neutral hydrogen, i.e., our 
signal. The measured power spectrum of the 21cm intensity P(k) is then, 

< ¥g >V»=P(k)6fg =plP HI (k)S^ + g 

P(k)=p 2 P HI (k) + V R <J 2 pN , (A7) 

where Vr is the volume of a pixel. 
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